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We study the effects of phonons in iV-leg Hubbard ladders within the framework of a one-loop 
renormalization group. In particular, we explicitly demonstrate that the role of phonons changes 
qualitatively even in the simplest two-leg ladder, as compared to the single-chain system where 
phonons always dominate. Our numerical results suggest that in the spin-gapped phase of the 
two-leg ladder, the opening of the spin gap by electron-electron interaction also drives the electron- 
phonon interaction to strong coupling, but in a subdominant fashion. Therefore, even though the 
inclusion of phonons does not alter the phase, their subdominant relevance strongly renormalizes 
some physical properties below the energy scale of the spin gap. This might shine some light on the 
recent experiments showing an anomalous isotope effects in high-temperature superconductors. 
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How relevant phonons are to the high-temperature su- 
perconductors has been an issue of considerable contro- 
versy. At one extreme, it is believed that phonons are 
irrelevant to the physics of the cuprates which is dic- 
tated by electronic correlations. At the other, phonons 
are thought to be the glue for the Cooper pairing. Recent 
experiments of angle-resolved photoemission experiments 
indicate the significance of electron-phonon coupling in 
the cuprates 0,12, El and spark the interests in the con- 
troversial role played by phonons again. This is particu- 
larly so in the isotope- ARPES experiment by Gweon et 
al.[3 There, the effect of isotope substitution is signifi- 
cantly enhanced in the superconducting state where the 
spin gap has opened. 

The purpose of this work is to shed some light on the 
effect of electron-phonon couplings in correlated electron 
systems. The challenge here is to do a calculation which 
is controlled when appropriate limits are taken, and at 
the same time is unbiased in the sense that no partic- 
ular scenario is favored or implicitly assumed, e.g., by 
a certain choice of mean field. For the particular case 
of the high-T c materials, no calculational method is cur- 
rently available that would satisfy both criteria. In the 
broader context of correlated electron systems, however, 
we feel that the class of ladder systems especially lends 
itself to an investigation of the issues stated above: Only 
in ladder systems it is known that purely repulsive elec- 
tronic interactions can lead to a spin-gapped phase (with 
approximate <i-wave symmetry), in resemblance to sce- 
narios that have been proposed for the superconducting 
cuprates. If anything, this choice of system will take the 
bias off the phonons as the driving mechanism behind 
such a phase, since it already exists without their help. 
Moreover, this nontrivial, strongly correlated phase with 
various competing order parameters^ can be accessed 
at arbitrarily weak strength of the electronic couplings, 



making it possible to study this problem by means of the 
perturbative renormalization group (RG). 

While the equations we derive are valid for general N- 
leg Hubbard ladders, we focus on the simplest two-leg 
ladder in this paper. Our results demonstrate a clear 
distinction of the interplay between electronic and lat- 
tice degrees of freedom, depending on whether a finite 
spin gap develops in the ground states. In the spin- 
gapped phase, the electron-phonon couplings are always 
relevant. In particular, they are driven to large values 
at the energy scale of the spin gap, which is set by the 
electronic interactions. This happens even in the limit 
of infinitesimal bare values of electron-phonon couplings. 
In gapless phases, on the other hand, electron-phonon 
couplings are predicted to be irrelevant at least for suffi- 
ciently small bare strength. This dichotomy indeed bears 
some resemblance to the experiments on the cuprates. Q 
At the same time, it is found that the electron-phonon 
couplings, although relevant, remain subdominant com- 
pared to interelectronic couplings. In simple terms, while 
the phonons do not drive the spin-gapped phase, they do 
dress it up. We note that this result is well distinct from 
that for the spin-gapped phase in the single chain, where 
even very weak attractive retarded couplings (resulting 
from electron-phonon interactions) will eventually domi- 
nate over the nonretarded (electronic) ones and play the 
essential role of generating the spin gap.[5j Our findings 
in the Hubbard ladder thus seem to set an interesting 
balance between the two extremes mentioned previously. 
On the one hand, phonons are irrelevant in regard to the 
pairing mechanism. On the other, they do strongly in- 
fluence what we observe in experiments, e.g., the isotope 
effects. While the above picture is suggestive for future 
explorations in the cuprates, we do warn the readers that 
our analysis is confined to ladder systems, which cannot 
be straightforwardly generalized to correlated systems in 
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higher dimensions. 

We now turn to the rigorous renormalization group 
analysis of ladder systems. In this Communication, we 
consider the effect of phonons in the iV-chain Hubbard 
ladder by including retarded interactions into the RG 
treatment. In the following we shall assume general in- 
commensurate band filling and weak interactions. Un- 
der these conditions, the effective field theory for the 
iV-chain ladder is a 2iV-species massless Dirac theory 
with four-fcrmion interactions. The Hamiltonian density 
H = H +Hi of the effective theory is given by HHSIJ: 

ia. 
ija/3 

where open boundary conditions transverse to the chain 



direction are assumed. Here, R and L stand for right and 
left moving electrons,respectively, i and j are the band 
indices, and a, (3 are the spin indices. The superscripts 
£ and s of the interaction parameters refer to large and 
small momentum transfer scatterings. Following Ref.jfJ, 
the letters / ("forward") and c ("Cooper") are used to 
distinguish whether the band indices are conserved at the 
scattering vertices. For i = j, the two become identical, 
and we choose c^f = to avoid double counting (dif- 
ferent from the choice in Ref|y). All nonzero interaction 
parameters are symmetric in band indices. In addition, 
terms involving only the R or the L fermion operators 
do not participate in the RG at one-loop order, hence 
they are omitted. For the Hubbard ladder, the bare val- 
ues of all interaction parameters are equal to the on-site 
repulsion U. The one-loop RG recursion relations for the 
interaction parameters in Eq. Q are known, 
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where Cij S — c if + ^ijfif f° r notational convenience. 
In the above di = d/dl, with I being the logarithm of 
the ratio between the bare momentum cutoff Ao and 
the running cutoff A. Dimensionless couplings = 
g~ij/ir(vi + Vj) have been introduced with g = f e , f s ,c e ,c s 
and a! li:fc = (vi + v k ){v.j + v k )/[2v k (vi + Vj)]. 

The numerical integration of Eq. (j5J shows that quite 
generally for the iV-chain system, |(j there exists a scale Id 
where the RG flows diverge due to the opening of spin or 
charge gaps in some/all partially filled bands. To ensure 
the applicability of the perturbative results, the flows are 
usually terminated at a scale l c (close to but less than 
Id), where the largest couplings become of order I. Phys- 
ically, the scale l c is identified with the energy scale at 
which charge/spin gaps open. 

To account for the effects of phonons we include re- 
tarded four- fermion interactions in Eq. The latter 
are effective interactions which are formally identical to 
those displayed in Eq. (Q), except that a cutoff wd is 
imposed for frequency transfer at the retarded vertices. 
By standard momentum shell RG,Ei| one can derive RG 
equations for both nonretarded and retarded couplings 
before the momentum cutoff A reaches ujp ly at a scale 
l = l D , using the method of Refs. El El El E3 Below 
this scale, the retarded couplings are then regarded as 



additive renormalizations to the nonretarded ones. 

For I < l D , the following flow equations are obtained 
for the retarded interactions: 
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whereas the nonretarded flow equations J2J) remain un- 
altered. In the above, we use the notation g for the re- 
tarded version of the dimensionless coupling constant g. 
Note that the retarded couplings with small momentum 
transfer /A , gy do not flow at one-loop order. 

We will assume in the following that Eqs. 10) remain 
valid until the energy scale of the electronic spin/charge 
gaps is reached. Physically, this corresponds to the as- 
sumption that the Debye scale lo d is of the order of the 
electronic gaps or less. The evolution of the interactions 
will then be governed by the coupled set 10) and (0) 
of nonlinear differential equations as higher-energy elec- 
tronic excitations are recursively integrated out. Since 
the nonretarded flow equations (0) are unaffected by 
the retarded couplings, we may regard their solutions as 
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known. Thus, all we need to do is to substitute them 
into Eqs. J2J and solve for the renormalized retarded 
interaction parameters. 

It is rather remarkable that the RG equations © for 
retarded couplings are completey decoupled by introduc- 
ing the following linear combinations of the original cou- 
plings, hfj = ffj±cfj. After the change of coupling basis, 
the RG equations take the simple decoupled form, 



dl~hti=pt j hf j -2(h±) 



where p i 



(2/f.-4/£)±(2cf,-4cfA 



(4) 

Note that, since 
Cu = 0, the label ± in Eq. Q is redundant for i — j and 
thus will be dropped. Omitting indices, the solution of 
Eq. igj is of the form 



h(l) = M(l) / dl'2M{l') + ^ 
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where M(l) — exp (j* dl' . In the vicinity of the 

divergent scale Id where the nonretarded interactions di- 
verge, the coefficient function p(l) diverges as 



P(l) 



(6) 



where the coefficient [3 depends on the subscript ij and 
the superscript ± of p. Because of the above divergence 
the retarded couplings will be driven to large values. 
Given Eq. ©, the following three cases need to be dis- 
tinguished: 

(1) (3 > 1. The integral over M in Eq. 10 diverges. 
The retarded couplings will then diverge as (Id — 1 
with Id < Id- In the interesting case where the retarded 
interaction is attractive h(0) < (which is realized for 
hfj and ha), the quantity in the parenthesis in Eq. (J5J 

vanishes before the scale Id is reached, and hence Id < 
Id- In this case, retarded couplings eventually dominate 
over nonretarded ones, so long as uj d is sufficiently small 
(compared with gaps in the system). For h(0) > we 
have ld = h- 

(2) < (3 < 1. A finite critical initial coupling 



dl2M(l) 



< 



(7) 



can be identified. For h(0) < h c < 0, the behavior is sim- 
ilar to case 1) with Id < Id- For h(0) > h c , h(l) diverges 
as (Id — . This implies that h is still relevant, but 
is ultimately subdominant compared to the nonretarded 
couplings, which diverge as (Id — ■ 

(3) (3 < 0. Again, the behavior for h(Q) < h c is as in 
case 1) with Id < Id, but otherwise h(l) is irrelevant and 
vanishes as (I — Z^)'' 3 '- 

At each point of the CI SO phase of the pure Hubbard 
ladder, the retarded couplings may now be classified ac- 
cording to the above scheme. We shall now do so for the 
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FIG. 1: (3 parameters defined in Eq. for the retarded 
couplings of a two-leg ladder. Note that the intraband cou- 
pling hn divides the C1S0 phase into two regimes regarding 
whether phonons play the dominant role in determining the 
gaps. Numerically, j3 has been obtained at l c for U/t = 10 -5 . 
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FIG. 2: Phase diagram of the two-leg Hubbard ladder as a 
function of t±/t and particle density. t± and t are hopping 
amplitudes along rungs and chains. Phases are denoted by 
CmSn, where m and n are the number of gapless charge and 
spin modes, respectively. In the regime below the thick line, 
both bands are partially occupied. The dash-dotted line in- 
dicates where umklapp processes are important (Refs. 0, 0) 
which were neglected in this work. The velocity ratio V1/V2 
equals 1 on the t± = and n = 1 axes, and diverges towards 
the C1S1 phase boundary. 



two-leg ladder. The presentation of the results is greatly 
simplified by the fact that the parameter (3 in Eq. © only 
depends on the ratio of the two Fermi velocities V1/V2, 
where the indices 1 and 2 refer to the majority (bonding) 
and minority (antibonding) bands, respectively. 

The (3 parameters of the retarded couplings according 
to Eq. © are displayed in Fig. ^ As is readily seen, 
below v\/v2 — 6 in the C1S0 phase, the (3 values for the 
most relevant retarded couplings behave according to the 
second case discussed above. The critical initial values, 
h c , for these couplings defined in Eq. (J7J are of the order 
of the bare electronic on-site repulsion J7.^(| I n mate- 
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rials where the bare retarded couplings are small com- 
pared to U, the majority of the fully spin-gapped C1S0 
phase (Fig. |2J will then have retarded couplings diverge 
subdominantly, implying that they will not surpass their 
electronic counterparts. This finding is dramatically dif- 
ferent from the conventional wisdom drawn from strictly 
one-dimensional systems. |l2t 1 1 3| where even very weak 
retarded couplings will dominate the nonretarded ones 
and dictate the phase. On the other hand, the effects of 
phonons are, though less dominant, by no means negligi- 
ble in the present scenario. This is particularly true for 
physical effects, such as the isotope effect, that primarily 
depend on the strength of the electron-phonon interac- 
tion: The renormalization of electronic properties due 
to phonons will be considerably stronger under circum- 
stances where the electronic interactions are relevant and 
open up gaps in the spectrum. This is best demonstrated 
by comparing the behavior described above for the CI SO 
phase to that found in the C2S2 phase, where electronic 
interactions are irrelevant. Note that, in the C2S2 phase 
(no spin gap), the quantity p(l) will not show the behav- 
ior displayed in Eq. J^J, but rather approach a constant 
value as I — ^ oo. In this case, even though the retarded 
couplings may formally diverge for p(oo) > 0, the energy 
scale associated with this divergence vanishes as a power 
of the bare electron-phonon couplings. 16] They are thus 
irrelevant when the cutoff lo d is taken into account, below 
which Eqs. @ are not valid. Again we note that this di- 
chotomy between gapped and gapless phases is paralleled 
by the experiment of Ref. 0, where the reported isotope 
effect is much more pronounced for temperatures and in 
regions of k space where a d-wave superconducting gap 
prevails. 

The presence of phonons also has an influence on the 
spin gaps. First, we point out that in simple models the 
bare value of is usually zero owing to time-reversal 
symmetry. From Eq. fll)l. h^ 2 will then remain zero 
within one-loop RG. hj 2 is irrelevant for small initial 
values due to the negative sign of the corresponding 
parameter (Fig. Under these conditions the relevant 
retarded couplings are intraband couplings. The latter 
become large and negative under RG, whereas their elec- 
tronic, initially repulsive, counterparts change sign un- 
der RG and eventually also flow to large negative values. 
Hence attractive retarded interactions tend to enhance 
the effects of electronic interaction in causing the spin 
gap. In addition, the increase of attractive interactions 
at large momentum transfer due to electron-phonon cou- 
pling will tend to enhance charge-density wave (CDW) 
correlations, similar to the case of the single-chain sys- 
tem. 

In a slim regime of the phase diagram (above v i /v 2 ~ 
6), the fi value of h\\ is larger than 1, such that the first 
case mentioned above applies. In this narrow regime, 
the phonons take over the lead as in the single chain. 
This change can be roughly understood as the follows: 



Due to the large mismatch in Fermi velocities, the two 
bands decouple as two independent single-chain systems. 
In fact, in the regime above v\jvi ~ 8 (often misidenti- 
fied as C2S1 phase in the literature), fin « 3/2, which is 
the same as in the spin-gapped phase of a single chain. 
It is a bit tricky to determine the phase in this regime. 
Naively, at the cutoff scale I — l c , the dominant retarded 
coupling hn(l c ) reaches order 1 and opens up the spin 
gap in the bonding band, while all other retarded and 
nonretarded couplings remain small. It is tempting to 
reach the wrong conclusion that the phase C2S1 is re- 
vived. Upon more careful analysis, 0, the subdomi- 
nant couplings ft.22 , /I2 & lso develop a minor spin gap in 
the antibonding band. In addition, the subdominant cou- 
pling c[ 2 > pins the relative phases between two bands 
and leads to the "<i-wave" pairing symmetry. Therefore, 
the final phase is again the same C1S0, except the spin 
gap in the bonding band is much larger than the others. 

To conclude, we have developed a general RG scheme 
to analyze the effects of phonons in the correlated iV-leg 
ladder. In particular, we demonstrate that the phonon- 
mediated attractive interaction, although divergent, re- 
mains subdominant compared to the electron-electron in- 
teractions in the two-leg ladder. This result is particu- 
larly interesting when compared to that for the Luther- 
Emery phase of a single chain, where phonon effects can 
be even more profound. The general tendency suggested 
by these findings is that as one moves away from a purely 
one-dimensional chain, the essential part in determining 
the phase diagram is played by instantaneous interelec- 
tronic interactions. Phonons may, however, inject an 
important bias when two leading competing electronic 
divergences are equal in strength. At the same time, 
the results obtained here highlight that electron-phonon 
couplings may have an enhanced influence on quantita- 
tive aspects, in particular for measurements that probe 
energies below the spin gap set by electronic interactions. 
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